Crystal water induced switching of magnetically active orbitals in CuCl2 
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The dehydration of CuCl2-2H20 to CuCl2 leads to a dramatic change in magnetic behavior 
and ground state. Combining density functional electronic structure and model calculations with 
thermodynamical measurements we reveal the microscopic origin of this unexpected incident - 
a crystal water driven switching of the magnetically active orbitals. This switching results in a 
fundamental change of the coupling regime from a three-dimensional antiferromagnet to a quasi 
one-dimensional behavior. CuCl2 can be well described as a frustrated J1-J2 Ifeisenberg chain with 
ferromagnetic exchange Ji and J2/J1 ~ —1.5 for which a helical ground state is predicted. 

PACS numbers: 



I. INTRODUCTION 

Low-dimensional spin 1/2 magnets are of wide interest 
in solid state physics since they are ideal objects to study 
the interplay of dimensionality, magnetic frustration and 
strong quantum fluctuations. These compounds can be 
described often very successfully based on their magnet- 
ically active structural building blocks and their linking. 
Typical examples for such building blocks are Ti(III)06 
octahedra, V(IV)05 square pyramids or Cu(II)04 pla- 
quettes that form various spin 1/2 networks like quasi 
one-dimensional (ID) chains or ladders. Nevertheless, 
for a reliable and accurate description of such networks, 
precise model parameters are a precondition, especially 
in the vicinity of quantum critical points. However, since 
for new compounds these parameters are unknown, it is 
common to transfer the known parameters from related, 
similar systems in a slightly renormalized form according 
to changed distances and/or bond angles. 

A typical example for the application of this strat- 
egy are compounds that contain crystal water in differ- 
ent amounts. Although in some cases the topology of 
the magnetic network and the related magnetic proper- 
ties totally change upon dehydration like in the case of 
CuSi03-H2C>i''^ and CuSiOa'^''', for most compounds only 
moderate structural changes with respect to the magnetic 
network are observed. It is generally assumed that in 
this case crystal water leads mainly to a modest change 
of the crystal field for the magnetic ion. In turn, small 
changes in the crystal field only, would directly suggest a 
description within the same model with slightly revised 
parametersi^ This leads to the common believe that crys- 
tal water plays only a minor role regarding the magnetic 
properties for compounds where crystal structure is ba- 
sically preserved upon water intercalation. 

Here, in contrast, we show that the hydration of 
CuCl2 to CuCl2-2H20 fundamentally changes the mag- 
netic properties, although the topology of the cova- 
lent Cu-Cl network is seemingly unchanged. Whereas 
CuCl2-2H20 is a classical three-dimensional (3D) anti- 
ferromagnet (AFM) with a Neel temperature of 4.3 K, 
we establish the dehydrated species as an example for 



a quasi ID chain compound. The results of susceptibil- 
ity measurements, density functional and model calcula- 
tions can be consistently understood from a reorientation 
of the magnetically active Cu orbital driven by the hy- 
dration. Regarding the well known 3D magnetic nature 
of CuCl2-2H20f^'i the quasi ID behavior of the water 
free compound CuCl2 is rather surprising: CuCl2 is a 
J1-J2 Heisenberg chain model compound with FM near- 
est neighbor (NN) exchange Ji and AFM next-nearest 
neighbor (NNN) exchange J2- According to the esti- 
mated J2/J1 ^ —1.5 we predict a helical magnetic order 
below the observed transition at 24 K. Earlier studies^.1^ 
that tried to model CuCl2 as a spin 1/2 chain found con- 
siderable deviations from a ID behavior since in these 
investigations only NN AFM coupling was considered. 
Alternatively a ID model with AFM NN and NNN ex- 
changes was suggested, but this model yields far too large 
couplings.— 



II. CRYSTAL STRUCTURE 

The crystal structure of CuCl2 is presented in Fig. [TJ 
Cu and CI form a covalent network of edge shared CuCU 
plaquettes running along b direction. Such fourfold pla- 
nar Cu^"*" coordination suggests a strong analogy to the 
undoped cuprates. The Cu-Cl-Cu bond angle is 93.6° 
and thus very similar to that in the Cu02 chain cuprate 
family, where bond angles close to 90° result in FM 
NN exchange like in Li2Cu02,"'^2 LiCusO^i^iiiii^ and 
Li2ZrCu042^J^ As in the latter two compounds, these 
chains are arranged in layers (Fig. [Tl top), suggesting a 
rather weak exchange between the layers. Since even the 
arrangement of the chains within the layers is very simi- 
lar, a quasi ID behavior might be expected from a mere 
comparison with these cuprate crystal structures. 

When CuCl2 is exposed to moisture, H2O enters the 
space between the chain layers, finally forming the fully 
hydrated CuCl2-2H20 (Fig. [3 bottom). Although the 
crystal structure seems very similar at a first glance, the 
crystal water induces several changes: (i) The inter-layer 
distance increases, (ii) the CuCl2 chains shift with re- 
spect to each other and, (iii) the Cu-Cl distances within 
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FIG. 1: (color online) Crystal structures of CuCl2 (top) and 
CuCb -21120 (bottom). The shortest Cu-Cl bonds are high- 
lighted (brown) . The crystal water is intercalated between the 
CuCl2 layers. The dashed line in the CuCl2-2H20 structure 
marks the unit cell of the water- free CuCl2. 



the chains are modified. Whereas the structural changes 
can be easily understood by packing and electrostatics 
(negatively polarized is situated close to the Cu^"*" 
ions, H+ is attracted by Cl~), the origin of the dras- 
tic change of magnetic properties - ID versus 3D - is 
far from obvious. Unraveling the underlying microscopic 
physics is the aim of our joint theoretical and experimen- 
tal study. 



III. METHODS 

Polycrystalline CuCl2 was prepared by dehydration 
of CuCl2-2H20 (Alfa Aesar 99.999%) under vacuum at 
390 K. Single crystals were grown by chemical trans- 
port in a temperature gradient from 650 K to 575 K with 
AICI3 (Alfa Aesar 99.999%, ultra dry) as transport agent. 
The chemical characterization of CuCl2 -21120 and of the 
CuCl2 crystals was carried out by X-ray powder diffrac- 
tion, DSC/TG-niethods and chemical analysis. The heat 
of dehydration was determined by DSCM Based on five 
independent measurements, the heat of dehydration is 
^^dchyd. = (117 ± 2)kJ/mol at 400 K. 

Magnetization was measured in a SQUID magnetome- 
ter (1.8 — 300 K) in magnetic fields up to 7T. Heat capac- 
ity was determined by a relaxation method in the same 
temperature range up to ^qH = 9T. 

Exact diagonalization (ED) of the J1-J2 Heisenberg 
Hamiltonian has been performed on TV = 16 sites clus- 
ters using the ALPS code.l^ The low-temperature behav- 
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FIG. 2: (color online) Magnetic susceptibility of single crys- 
talline CuCl2 (top) and CUCI2-2H2O powder (bottom) as a 
function of temperature {fioH = IT). TMRG fits are given 
by the dashed lines. The insets show the derivative dx/dT, 
the ordering temperatures are indicated by sharp peaks. 
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FIG. 3: Temperature dependence of the specific heat Cp/T 
of single crystalline CuCl2 . 



ior of the magnetic susceptibility has been simulated us- 
ing the transfer-matrix density-matrix renormalization- 
group (TMRG) method^ 

For the electronic structure calculations the full- 
potential local-orbital scheme FPLO (version: fploT.OO- 
28) within the local (spin) density approximation 
(L(S)DA) was used.— In the scalar relativistic calcula- 
tions the exchange and correlation potential of Perdew 
and Wang was chosen.^- To consider the strong elec- 
tron correlations for the Cu configuration, we use 
the LSDA-|-J7^^ approximation varying Ud in the phys- 
ically relevant range from 6 ~ 8.5 eV. The LDA results 
were mapped onto an effective tight-binding model (TB) 
and subsequently to a Hubbard and a Heisenberg model. 
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IV. RESULTS AND DISCUSSION 
A. Thermodynamic measurements 

The susceptibility data for both compounds are shown 
in Fig. [21 CuCl2 exhibits a broad maximum at T,„ax ~ 
75 K as a fingerprint of quasi fD behavior—. An AFM 
Curie- Weiss temperature Qqw = +107 K has been ex- 
tracted from the high temperature region. A sharp kink 
at Tjv = 24K (see Fig. [2] and inset dx/dT) followed by a 
rapid drop of x indicates a magnetic phase transition as 
earlier suggested.^ 

The measured zero-field specific heat as a function of 
temperature for CuCl2 is shown in Fig. [3l Our data 
agree well with earlier studiesi^iiS The specific heat curve 
shows a pronounced lambda shape anomaly at = 24 K 
due to the onset of long range AFM ordering. The 
ordering temperature from Cp is in perfect agreement 
with Tjv evaluated from susceptibility. Well below T/v 
the total specific heat is described by Cp = (3T^ with 
13 = 0.8514(3)mJ mor^ K^i (fit for T < lOK), indicat- 
ing that both phononic and magnetic contributions to 
Cp{T) are oc T^. The data in a field hqH = 9T show no 
visible differences to the zero-field data. 

In contrast to the quasi ID susceptibility of CuCl2, 
the hydrated system shows an increasing x down to low 
temperatures right above the AFM phase transition at 
4.3 K (see Fig. [2] and inset dx/dT) in perfect agreement 
with earlier measurements A Curie- Weiss temperature 
©cw — -t-5.3K indicating weak AFM interactions has 
been evaluated. 



B. Band structure calculations 

The essentially different character of the susceptibil- 
ity of both compounds points to a changed coupling 
regime rather than to a mere re-scaling according to the 
modified atomic distances. To construct an appropri- 
ate microscopic model based on the relevant interactions 
we perform ab-initio electronic structure calculations, as 
successfully demonstrated earlier for the closely related 
Cu02 chain compound famil y^^'^^'^^i^^ . 

Total and partial densities of states (DOS) for both 
compounds are pictured in Fig. H) On a coarse energy 
scale both systems are similar, the contribution of the ad- 
ditional H states to the valence region is negligible. Both 
compounds show half-filled, well separated anti-bonding 
bands at the Fermi level. This metallic behavior is in 
contrast to the experiment and is a well known short- 
coming of the LDA due to the underestimation of the 
strong Coulomb repulsion for the Cu^+ Sd^ configura- 
tion. The observed insulating ground state is obtained (i) 
within the LDA-|-?7 approximation or (ii) by a model ap- 
proach mapping the relevant low lying LDA states onto 
an effective TB model, and, including the correlations, 
subsequently onto a Hubbard and a Heisenberg model. 
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FIG. 4: (color online) Total and partial density of states and 
band structure of the anti-bonding band of CuCb (a,e) and 
CuCl2 -2112 (b,f). The insets show the orbital character of 
the anti-bonding Cu-Cl (a) and Cu-Cl-0 (b) states, respec- 
tively. The bandwidths are compared in (d). The fits (full 
lines (e,f)) of the TB model (c) are superimposed to the LDA 
band structures (circles (e,f)). The two bands in (f) originate 
from the two Cu atoms per cell in CuCl2-2H20. 



A closer inspection of the DOS and the related band 
structure reveals two important differences: (i) Whereas 
the width of the anti-bonding band in CuCl2 is 0.8 eV - 
rather typical for ID edge-shared Cu02 chains'^^ - the 
bandwidth in CuCl2 -21120 is reduced by more than a 
factor of three to about 0.25 eV (Fig. [JJi). (ii) The mag- 
netically active anti-bonding band in CuCl2 is formed 
exclusively by Cu-Cl dpa states (inset Fig. |4^) corre- 
sponding to the bonds pictured in Fig. [TJ In contrast, 
for CuCl2 -2112 the O 2p orbitals that are directed to- 
wards the Cu contribute significantly to this band (inset 
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FIG. 5: (color online) Sketch of the magnetically active 
pdu-plaquette for CuCl2 (top) and CuCl2-2H20 (bottom). 
Whereas these orbitals form quasi ID edge shared chains in 
CuCb, the plaquettes in CuCb -21120 are isolated resulting 
in weak, but 3D interactions. 

Fig. [Dd) . This leads to the formation of a new dpa orbital 
perpendicular to the original ones. The related CuCU 
or CUCI2O2 plaquettes are shown in Fig. O These pla- 
quettes, relevant for the magnetic couplings, form edge- 
shared chains in CuCl2, whereas they are isolated in the 
hydrated system. 

Naturally, this "orbital switching" induced by the crys- 
tal water implies a change of the coupling regime from ID 
to 3D. To study these changes on a quantitative level, we 
constructed a TB model (Fig. Uj:) for both compounds 
and fitted it to the relevant LDA bands (Fig. |l^,f). 
The leading transfer integrals ti are given in Table [J 
The calculated U confirm the intuitive picture that cor- 
responds to Fig. El CuCl2 shows quasi ID dispersion 
along the chain with dominating NNN hopping t2, and a 
considerably weaker coupling between the chains, while 
the coupling between adjacent layers is very small. On 
the other hand, the changed plaquette arrangement in 
CuCl2 -21120, induced by the crystal water, leads to a 
strongly reduced band width (Fig. |4]i) and correspond- 
ingly small isotropic (3D) transfer integrals (Table 

For the strongly correlated limit (C/off S> ti) at half fill- 
ing, the TB model can be mapped via a Hubbard model 
to a Heisenberg model with resulting AFM exchange cou- 
plings Ji = Aif/UeS where UeS is the correlation in the 
effective one-band description. Depending on the choice 
of UeS within a reasonable rang o^i^^'^^ (C^cff = 3.5 — 4 
eV) this leads to exchange constants = 3 — 5K for 
the three leading couplings in CUCI2-2H2O (Table 
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TABLE I: Calculated leading hopping integrals ti for the ef- 
fective TB model shown in Fig. |4};. 



These J's are perfectly in line with the experimentally 
observed Tat = 4.3 K as could be expected for an almost 
isotropic 3D coupling. For CuCl2, the leading NNN t2 
results in J2 — 160 — 180 K. Although of the correct or- 
der of magnitude compared to Tmax and Ocw , this would 
clearly exceed the "overall AFM coupling" in the com- 
pound without additional FM interactions. Sizable FM 
interactions are typical for close to 90° bond angles ac- 
cording to the Goodenough-Kanamori- Anderson rules^^ 
and therefore expected for the NN J\ in CuCl2. For 
a quantitative estimate of the FM contributions to the 
leading J's we apply LSDA-ffJ calculations for different 
spin arrangements in magnetic super cells. 

Mapping the resulting total energy differences to the 
Heisenberg model, we obtain the following total exchange 
integrals (C/d = 7 ± O.SeV)^^. ^ -(150±10)K, 
J2 = 155 ± 25K and J'f = 35 ± 5K. Within the error 
bars, the latter two agree very well with the J values cal- 
culated from the corresponding t's of the TB approach;^ 
indicating that FM contributions beyond NN are rather 
small. In contrast, we find a large FM contribution of 
about 175 K to Ji as expected from the Cu-Cl-Cu bond 
angle of 93.6°. The size of the FM contribution to J\ fits 
well to related edge-shared Cu02 chain compounds li^ii^ 

The resulting leading exchange interactions confirm 
the intuitive picture of a quasi ID chain model com- 
pound with small inter-chain coupling, very similar to 
LiCu202»^i^i^ Therefore, the magnetic ground state is 
mainly determined by the ratio a=J2/Ji of the frustrat- 
ing main interactions along the chains. For CuCl2, we 
find a = —(1.0 ± 0.1) and predict a ground state well in 
the helical ordered region of the Ji-J2-phase diagram. 

Thus, the dehydration of CUCI2-2H2O to CuCl2 leads 
to a drastic change of the coupling regime from 3D to 
quasi ID and a completely different ground state that 
can be traced back to a switch of the magnetically active 
orbital. 



C. Model analysis of magnetic susceptibility 

For an independent evaluation of the leading exchange 
interactions in CuCl2 we simulated x(T) within a spin- 
1/2 J1-J2 Heisenberg model for various J2/J1 ratios 
using the TMRG technique and fitted the resulting 
X*(T/J) curves to the measured x(^)- Rather typical 
for the J1-J2 Heisenberg model,— we find two possi- 
ble solutions for the fit: (i) a = -h3.0 with Ji = 120 K 
and J2 = 40 K (AFM solution) and (ii) a = -1.5 with 
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Ji = -90 K and J2 = 135 K (FM solution). The FM so- 
lution (ii) is in rather good agreement with the estimates 
from our ab-initio calculations. The corresponding fits 
are shown in Fig. [51 The AFM solution (i) can be dis- 
carded regarding our calculational results and the close 
to 90° Cu-O-Cu bond angle. 

In a naive approach, using the relation &cw ~ 
l/2(Ji + J2 + Jic),'^'^' the theoretically estimated 
Qcw ~ '^'^^ ±10 K seems to be inconsistent with the ex- 
perimental 8cw ~ -1-100 K. Thus, we choose a more so- 
phisticated procedure performing ED studies for the FM 
solution^ which yield the expected 9cw ~ l/2( Ji + J2), 
but only at high temperatures T > 10 J2. Since this 
high temperature region is inaccessible to experiments, 
we choose a temperature window of the Curie- Weiss 
fit in order to obtain Qcw for the highest tempera- 
tures measured (225 K < T < 300 K corresponding to 
I.7J2 <T < 2.2J2). Using Ji and J2 from the TMRG, 
we obtain Q*^^ = +72K. The remaining difference be- 
tween and of about 30 K corresponds in very 
good agreement to the inter-chain coupling Jjc of about 
35 K neglected in the ED simulations. 

D. Thermochemical properties 

The dramatic effect of the dehydration of CuCl2 -21120 
certainly raises the question whether this is really a typi- 
cal dehydration process or if the involvement of the crys- 
tal water related oxygen in the magnetic exchange via 
covalent Cu-0 bonds is an indication for a chemical re- 
action on a different energy scale. The dehydration en- 
thalpy AH^^^y^ = 117kJ/mol of CuCl2-2H20 is a typ- 
ical (small) value for this class of materialsi ^^'^^ Our 
ab-initio estimate for the dehydration enthalpy yields 
^^dchyd ~ 95kJ/mol which is in quite good agreement 
with the measured value, providing additional confidence 
to the reliability of the calculational procedure. 



V. SUMMARY 



In summary, our joint theoretical and experimental 
study provides a consistent explanation for the funda- 
mentally different magnetic properties of CuCl2 -21120 
and CuCl2. Whereas CUCI2-2H2O is a quite isotropic 3D 
AFM with small exchange couplings due to the orienta- 
tion of neighboring isolated CUCI2O2 plaquettes, CuCl2 
can be well understood in terms of a ID FM-AFM J1-J2 
chain model. This extension of the originally, critically 
discussed ID AFM-NN only modelS. re-establishes the 
pronounced ID nature of the magnetism in CUCI2. The 
dramatic change of magnetic properties between both 
compounds can be traced back to a switch of the mag- 
netically active orbital induced by crystal water. From 
our ab-initio calculations and model studies of the mea- 
sured susceptibility using the transfer-matrix density- 
matrix renormalization-group and exact diagonalization 
techniques we predict a helical ground state and likely re- 
lated multiferroic behavior for CuCl2 driven by strong in- 
chain frustration originating from FM nearest neighbor 
and AFM next-nearest neighbor exchange interactions. 
Our study reveals that crystal water can have crucial in- 
fluence on the electronic and magnetic properties of low 
dimensional magnets. More general, our work empha- 
sis that a transfer of model parameters from seemingly 
closely related systems is rather dangerous. The unaware 
neglect of this fact will lead to an at best inaccurate de- 
scription of the physical properties in many cases. 

We acknowledge R. Kremer for valuable discussions, 
T. Xiang for providing the TMRG code, and S. Miiller 
for supporting DTA and DSC studies. 
Note added in revision: While revising our manuscript, 
we learned that our magnetic model and predicted ground 
state was confirmed by an independent study^ 
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These compounds are known as J1-J2 chain models with 
FM NN and AFM NNN exchange, where the different 
J2 / Ji ratios lead to different ground states. 
Netzsch DSC 204 thermocouple, platinum crucible, heat- 
ing rate lOK/min, dry argon atmosphere, 10 mg sample. 
The absence of an impurity related Curie tail at low tem- 
peratures indicates the high quality of the sample. 
For comparison, CuGeOa shows a bandwidth of 0.95 eV 
for the anti-bonding band, CuSiOa about 0.65 eV and 
LiCuV04 about 0.7 eV. 

The chosen region for Ud covers the experimental NN ex- 
change Ji in La2Cu04. 

This good agreement justifies a posteriori the transfer of 
the Ueff and Ud values from the CUO2 chain compounds. 
Even for medium size clusters (A'^ = 16), ED perfectly 
describes the high-temperature region of x{T) which obeys 
a Curie- Weiss law. 
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